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1 Introduction 

Supersymmetry, an elegant extension of spacetime symmetry, is the leading candidate for 
the new physics unfolded in the Large Hadron Collider (LHC) at CERN. The two detectors 
ATLAS and CMS of the LHC have been collecting data at a much faster than expected and 
their recent data significantly extend the exclusion limits for super symmetric particles. But 
the latest such data have so far been interpreted by the experiment in only two different 
supersymmetry breaking models: the constrained minimal supersymmetric standard model 
(CMSSM) and a simplified model with only squarks and gluions and massless neutralinos. 

Other supersymmetry breaking models should be extensively analyzed in the era of 
the LHC. One of those to be analyzed is gauge mediated supersymmetry breaking (GMSB) 
model [1-9] where messenger fields, charged under the Standard Model gauge symmetry, 
mediate the breakdown of supersymmetry in the hidden sector to the MSSM sector. The 
soft masses in the visible sector arise from quantum effects of the messengers so the super- 
symmetry breaking scale of the visible sector is much lowered than the grand unified theory 
(GUT) scale compared with the gravity mediated supersymmetry breaking scenario. 
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As the LHC continues to collect experimental data, the precise studies on the physical 
parameters of the SUSY particles will become important. This requires one to do quantum 
loop calculations on the SUSY parameters. For instance, the gluino pole mass of the 
CMSSM has been considered up to two-loop order in ref. [10-12] while the neutralino and 
chargino pole masses in ref. [13, 14]. 

Our interests lie in the gluino pole mass of the GMSB model at two-loop order. Among 
the various different GMSB models we choose the Minimal Gauge Mediation (MGM) where 
a pair of messenger fields, fundamental and antifundamental under the SU{3)c gauge 
symmetry, mediate the supersymmetry breaking to the MSSM sector. The gluino pole 
mass of the GMSB model at two-loop order was first discussed in ref. [15] where the 
authors made a prediction that the NLO correction to the gluino pole mass is up to 10% 
of the LO pole mass. On the other hand, our prediction is, as shown later, 20% or even 
more. The difference between their and our predictions arises from how to handle the IR 
behavior of SUSY QCD. We state that our treatment of the IR behavior is more consistent 
with perturbative theory rather than theirs. We will rigorously discuss it in Section 5. 

For the renormalization of the MGM lagrangian parameters the DR scheme is adopted [16- 
18]. It is based on regularization by dimensional reduction along with modified minimal 
subtraction(MS) scheme. Not only the messenger fields wavefunctions but also their masses 
are renormalized. The MSSM quarks and squarks contribute to the gluino pole mass at 
two-loop order through the renormalization of the gluon(gluino) wavefunctions and the 
gauge coupling at one-loop order. We will also take these contributions into account. 

In this paper, we follow the two-component formalism to derive the self-energy func- 
tions in ref. [12, 19], and then present the analytic results of the self-energy functions up 
to two-loop order relevant to the gluino pole mass. We also perform a numerical analysis 
for the NLO correction of the gluino pole mass. 

2 Self-energy functions and pole masses for two-component spinors 

We briefiy review the self-energy functions for fermions in two-component notation and 
then describe how to compute the loop-corrected gluino pole mass of the MGM. All the 
details can be found in ref. [19] . 

We consider a theory with left-handed fermion degrees of freedom ipj with an index j = 
1, 2, • • • , A^. The self-energy functions for fermions in two-component notation are depicted 
in figure 1, where the shaded circles denote the sum of all one-particle irreducible(lPI), 
connected Feynman diagrams, and the external legs are amputated ^. In figure 1 the self- 
energy functions for two-component fermions are denoted by 3, H"^, fi, and 12, and the 
four-momentum p fiows from right to left. 

The pole mass is defined by the position of the complex pole in the propagator and is 
a gauge-invariant and renormalization scale-invariant quantity. The pole mass of a fermion 
can be found by considering its rest frame, in which the space components of the external 
momentum p^^ vanish. This reduces the spinor index dependence to a triviality. Setting 

^At this stage, counterterm corrections are not reckoned with. 
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Figure 1. The self-energy functions for two-component fermions. 

= 0), we search for the values of s. In other words, the poles of the full propagator 
(which are in general complex) 

Spoiej = {Mj - iTj/2f (2.1) 

are formally the solutions to the non-linear equation ^ 

det[sl - (1 - H^)~^(m-F n)(l - Sy^{rn + n)] = 0, (2.2) 

where is the symmetric N x N tree-level fermion mass matrix with mikfn^^ = mf6j. 

The gluino pole mass of the MGM is zero at tree-level and arises at one-loop order. 
Moreover the gluino do not mix other fermions so the master equation reduces to a simple 
equation as 

= 0. (2.3) 



(i-sy 

The solution for Eq. (2.3) is given as 

V~s = Y^, (2.4) 

and is perturbatively calculated as 

^ = n(') + [n(^)E('^ + + , (2.5) 

where the self-energy functions are expanded in powers of a^: 

J] = 0(^) + J](2) + ... ^ (2.6) 
E = + + . . . . (2.7) 

We use an iteration method to solve Eq. (2.5). We first get the leading order (LO) pole 
mass by substituting the tree level gluino mass (s = 0) to the one-loop function il(s)^"'^\ 
Then we substitute the LO pole mass into Eq. (2.5) to calculate the NLO pole mass. 
In order to evaluate the gluino pole mass up to two-loop order we need to evaluate the 
self-energy functions and E^^h 



Here Mj and Tj are a physical mass and width of a fermion. 
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3 Minimal gauge mediation 



3.1 Lagrangian of minimal gauge mediation 

The chiral superfields contained in the MGM are messengers ^ and ^, and a Goldstino 
multiplet X. The chiral superfield components corresponding to the chiral messenger field 
$ and <1>, in the fundamental and anti-fundamental representation of the SU{3)c gauge 
symmetry, are denoted as 

^ = {ct),i;,T), «I = ((A,^,-F). (3.1) 

The free lagrangian of the messenger fields and the SUSY Yang- Mills lagrangian are given 
as follows: 

- Vlgicb^T^il^ - i}T''4>^)X'' - v^gX^^Oj^T^cl) - c^T^V^t) 

(3.2) 

where (i) g is the gauge coupling, (ii) f"'^'^ are the antisymmetric structure constants of the 
gauge symmetry which satisfy 

[T'',T'']=if''''^T^, (3.3) 

for the generators T"" for the fundamental representation, (iii) A" is the gluino field, (iv) 
v^^^, the gluon field strength, 

<. = 5X-5-<-5r''«> (3-4) 
(v) T>^ is the real auxiliary boson field, and (vi) the covariant derivatives are defined as 

{D^~^y = d^~^'-igvl^{T'')^ 
{D,~^)\ = d^~^\ + igvl{T-)r^), 

(Z),,^), = a^V. + ^K(^")/^J' 
{D,,ij^\ = d^i,\ + igvl{T\^i^], 

D^X^ = d^X^ - gr^'^lX^. 
The Goldstino X couples to the messengers via a superpotential 

W = Xi^, (3.6) 

and has an expectation value: 

{X)=M + e'^Fx. (3.7) 

Its expectation value Fx sets the scale of SUSY breaking as \/Fx- The other expectation 
value M gives each messenger fermion a mass M, and the scalars mass squared masses 



(3.5) 
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Figure 2. The Feynman rules for the messenger fermions 

equal to it Fx- The corresponding mass eigenstates of the scalars are given in terms 
of the gauge eigenstates of the scalars as follows: 

For the Feynman rules we select the mass eigenstates for the messenger scalars. 
3.2 Feynman rules 

In order to systematically perform perturbative calculations we first establish a set of 
Feynman rules for the MGM using the two-component spinor formalism. Following the 
conventions for the Feynman rules in ref. [19] we acquire them as depicted in figure 2, 3, 4, 
and 5. Here we use the usual Feynman gauge for the gluon field. We omit the Feynman 
rules for the gluon for the sake of saving space. 

The messenger fermions have four different propagators: the two are chirality-preserving 
as shown in figure 2 (a) while the other two chirality-violating as shown in figure 2 (b) and 
(c). For the Feynman rules for chirality-preserving propagators we have two options for 
them: either a or a. For the Feynman rules for fermion-fermion-gluon vertices we can also 
select either or cr as shown in figure 2 (d) and (e). But the choice on a matrices for 
propagators and the vertices in a Feynman diagram should be simultaneously fixed. For 
instance, if one chooses a a for a propagator then its neighboring vertices must pick out a 
a. 

As for a gluino we have only a chirality-preserving propagator as shown in figure 3 
because it is massless. The gluino propagator as well as the gluino- gluino- gluon vertex can 
contain either a or cj like the messenger fermions. 

The messenger scalars have two different four-vertices which are obtained by integrat- 
ing out P^-term in Eq. (3.2): the same mass eigenstates have either the same directions of 
their arrows as shown in figure 4 (d) or the opposite directions of their arrows as shown in 
figure 4 (e). 
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Figure 3. The Feynman rules for the gluino 
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Figure 4. The Feynman rules for the messenger scalars 



There are eight different scalar-fermion- gluino three- vertices as shown in figure 5. One 
should pay attention to every direction of arrows of fields in the vertex diagrams. There 
are some important properties to recall: 

• A direction of arrow of a messenger fermion {ip or ip) is the same as that of gluino: 
either into a vertex or out of a vertex. 

• A direction of arrow of a messenger scalar is either the same with that of a messenger 
fermion tp or the opposite with that of the messenger fermion ip. 

• Among the eight vertices only the two in figure 5 (d) and (f) have a positive sign in 
the Feynman rules while the rest a negative sign. 

Although we do not include the MSSM quarks and squarks in the lagrangian (3.2) we 
need to take account of their effects on the renormalization procedure for later numerical 
analysis. Here we omit their Feynman rules which are referred to in ref. [19]. 
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Figure 5. The Feynman rules for the scalar-ferniion-gluino vertices. 



4 Renormalization 

We shall perform loop calculations for the self-energy and vertex functions using the Feyn- 
man rules as described in the previous section. We consider the one-loop and counterterm 
corrections to the propagators and vertices relevant for the gluino pole mass. As mentioned 
earlier we take into account the MSSM quarks and squarks whose propagators appear in 
the loops of figure 6, 7, and 13. 

All the corrections for the propagators are depicted in figure 6, 7, 8, 9, 10, and 11: 

• Figure 6 shows the one-loop and counterterm corrections to the gluino propagator. 
The solid lines represent either the messenger fermions or quarks while the dashed 
lines either the messenger scalars or squarks. The second to the last denotes the 
ghost loop corrections while the last stands for the counterterm. 

• Figure 7 shows the one-loop and counterterm corrections to the chirality-preserving 
gluino propagator. There are four different combinations of the messenger fermions 
and messenger scalars in the loop and their contributions are constructive, leading 
to UV divergence. We include not only the messenger fermion-scalar loops but also 
the quark-squark loops. 

• Figure 8 shows the one-loop corrections to the chirality- violating gluino propagator. 
There are four different combinations of the messenger fermions and messenger scalars 
in the loop but their contributions are destructive, resulting in no counterterm. They 
yield the pole mass of the gluino at one-loop order. 

• Figure 9 shows the one-loop and counterterm corrections to the chirality-preserving 
messenger fermion propagator. 
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• Figure 10 shows the one-loop and counterterm corrections to the chirality-violating 
messenger fermion propagator. 

• Figure 11 shows the one- loop and counterterm corrections to the messenger scalar 
propagator. 

All the corrections to the three-vertices are depicted in figure 12, 13, 14, and 15: 

• Figure 12 shows the one-loop and counterterm corrections to the messenger scalar- 
messenger scalar-gluon vertex. Reversing the arrow direction of either the messenger 
fermion or gluino is also taken into account. 

• Figure 13 shows the one-loop and counterterm corrections to the gluino-gluino-gluon 
vertex. We include not only the messenger fermion-scalar loops but also the quark- 
squark loops. 

• Figure 14 shows the one-loop and counterterm corrections to the messenger fermion- 
messenger fermion-gluon vertex. 

• Figure 15 shows the one-loop and counterterm corrections to the messenger fermion- 
messenger scalar-gluino vertex. 

The DR scheme treats the UV divergences in the same way as the MS scheme. The 
main difference between the two schemes is that degrees of freedom of a chiral fermion field 
in the DR scheme are set to 2 while those in the MS scheme to d/2. Therefore gauge fields 
in the DR scheme are accompanied by the e-scalar fields to maintain supersymmetry. At 
one-loop order, it is equivalent to set the dimensions in which the a matrices reside to 2. 

The renormalization Z factors are defined as 

and cancel off the inverse powers of e of the one-loop integrals in the figures. Both the 
one-loop corrections to the propagators and vertices yield the leading term of power series 
in Us for Z^^^'s. We list them in figure 16 and 17 where Nf is the number of MSSM quark 
flavor and A^mcss is the number of the messenger pairs. The Casimir operators C{R), C2{R) 
and C2{G) are defined as 

(T'^T'^)/ = C2(i?)<5/ (4.2) 
Tr[T'^T^] = C(i?)(5"^ (4.3) 

j^acdjbcd ^ C'2(G)(^"^ (4.4) 

whose values are C{R) = 1/2, C2{R) = 4/3 and C2{G) = 3 for the S'C/(3)c gauge symmetry, 
respectively. 

As a consistency check we evaluate Slavnov- Taylor identities, 

^ Zip^y Zfjy^y Z\\y Z^lj)l, Zj]J]J] 

Z^y/ Zy Zfpy/ Zy Z\-\J Zy ^ Z Zy ^ 



-8- 







'oooooootooooooooo^ 




'ooooooocStteooooooo 




ghost 

'cxxx»i jOTCcr 



Figure 6. The one-loop and counterterm corrections to tlie gluon propagator. 
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Figure 7. The one-loop and counterterm corrections to the chirality-preserving gluino propagator. 
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Figure 8. The one-loop corrections to the chirality-violating gluino propagator. They contributes 
to the pole mass of the gluino at one-loop order. 




Figure 9. The one-loop and counterterm corrections to the chirality-preserving messenger fermion 
propagator. 



where the various Z factors are calculated at one-loop order. We obtain the leading coef- 
ficient of the renormalization Z factor for the gauge coupling, 



= -^C2{G) + (Nf + iV,,ess)C(i?). 



(4.6) 
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Figure 10. The one-loop and counterterm corrections to the chirahty-violating messenger fermion 
propagator. 




Figure 11. The one- loop and counterterm corrections to the messenger scalar propagator. 




Figure 12. The one- loop and counterterm corrections to the messenger scalar-messenger scalar- 
gluon vertex. 

5 Self-energy functions 

5.1 Self-energy functions at one-loop order 

We first calculate the helicity-conserving self-energy function Q^^^ which arises from the 
one-loop diagrams as shown in figure 8. The sum of the four different configurations in the 
internal loops is free from UV and IR divergence, yielding the gluino pole mass at one-loop 
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Figure 13. The one-loop and countcrtcrm corrections to tlic gluino-gluino-gluon vertex. 




Figure 14. The one-loop and counterterm corrections to the messenger fermion-messenger fermion- 
gluon vertex. 




Figure 15. The one-loop and counterterm corrections to the messenger scalar- messenger fermion- 
gluon vertex. 

order as, 

^ = mg = n'^^\s) (5.1) 
= 2iV^essM„ess C{R) ( - (r?i - 1) ln(7?2 X + 1) - im - 1) ln(r?i x + 1) - (r ^ r=)) , 
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Figure 16. The countcrtcrm corrections to the propagators. 
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Figure 17. The counterterm corrections to the vertices. 
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(5.2) 

(5.3) 

(5.4) 
(5.5) 



and is the external four momentum. Mmcss is messenger fermion mass (which is denoted 
by M in section 3) and A is SUSY breaking scale of the visible sector which is equivalent to 
Fx/Miness- It is noted that the self-energy function 0^^^ has no explicit dependence on the 
renormalization scale, ^. Since the relevant range of s is much lower than the messenger 
scale, i.e. x < 10~^, it is a good approximation to take the limit x — )■ 0. Then O^^^ is given 
by [20, 21], 

r ln(r) 



J7(i)(s)U^o = (^)AiV,,ess2C(i?) 



{r-iy 



+ (r o f) 



(5.6) 
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Figure 18. The two-loop contributions to tiie cliirality-violating gaugino propagator. Each shaded 
circle denotes the one-loop and counterterm corrections either to the propagators or to the vertices. 
The dashed line can be either 0+ or while the solid line either ^ or ip. 



The helicity- violating self-energy function stemming from the diagrams in figure 7 
is given as 

'r + X — 1 



C{R) 



as 
All 



r — 1 r ln(r) 
+ 2 + + ^ + In 



(r/i - 1) ln(?72 X -I- 1) -I- (772 - 1) ln(?7i x + l) 
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mess 



+ {r + 2Nf[ In ( I - 2 
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(5.7) 



As mentioned earlier, it consists of the three parts: (i) the messenger loops, (ii) the quark- 
squark loops, and (iii) the gluino-gluon loops. In the limit x — )• 0, it is reduced to 



Att J 



+2A',(ln(^)-2))+ft(G)(l„(^)-2) 



(5.8) 



We will comment on the IR behavior in H^^^ (s) in the next subsection. 



5.2 Self-energy function at two-loop order 

In this section we describe how to calculate the self-energy function ^}^'^\s) in detail. 
All the Feynman diagrams relevant to fi^^^ are shown in figure 18, where each shaded 
circle denotes the one-loop and counterterm corrections either to the propagators or to the 
vertices in Section 4. The sum of the Feynman diagrams contains no UV divergence due 
to the inclusion of the counterterm corrections at one-loop level. 

We take a few steps to calculate the Feynman diagrams. We first decompose a Feynman 
integral with momentum tensors in numerators into integral forms of scalar products using 
the Gram determinant. Reducing the scalar Feynman integrals to the master integrals 
we use Laporta's algorithm [22], which systematically applies several reduction methods, 
such as Passarino-Veltman reduction [23], integration-by-part method [24, 25] and Lorentz 
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Figure 19. The master integrals. A dashed line represents a massless propagator while a solid 
line a massive propagator whose mass is explicitly noted. A dot on a solid(dashed) line denotes a 
double propagator. 

invariance method [26] . All the reduction procedures have been executed with our in-house 
Mathematica code. 

The sixteen master integrals are cataloged in figure 19. The dashed lines represent 
massless propagators while the solid lines massive propagators whose masses are explicitly 
noted. A dot on a solid(dashed) line denotes a double propagator. To evaluate the master 
integrals we mainly use the Mellin-Barnes representation which is described in ref. [27]. 
The analytic forms of the master integrals are fully described in appendix A. 

Performing all the procedures we acquire the self-energy function 0(^^(s). We checked 
that all the UV divergent poles exactly cancel out. The complete form of Q^'^\s) is so 
lengthy that it is attached in Appendix B. In the limit x — >• 0, the explicit formula of il^^^ 
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r — 1 



ln(r) + ln(f) 



+ (r -H- f) 



(5.9) 



The first three diagrams (from left to right) on the second row shown in figure 18 contribute 
to the term proportional to C^iG) in Eq. (5.9) while all the diagrams except for the one 
at the intersection of the second row and third column contribute to the term proportional 
to C2{R) in Eq. (5.9). The term proportional to C2{G) contains a term that represents 
IR behavior, i.e. ln(— s///^)^ since the gluon and gluino in the propagators are massless. 
On the other hand, the term proportional to C2{R) have an large logarithm, / M"^^^^) , 
since we take the renormalization scale /i as the physical gluino mass scale. 

As a side note we comment on the large logarithm, ln(/i2/M^ggg), in both H^^)(s) and 
Q.^'^\s). For A^'moss = 1 the term with / M"^^^^) in H(^)(s) gives less than 0.2 even 
though the Mmess goes up to 10^ TeV due to the small coupling constant as/(47r) = 0.008 
at /X ~ 1 TeV. Therefore we can safely do perturbation for the computation of the pole 
mass. However, this term increases linearly with A^mcss- For the case A^mess ^ 3 with large 
Mmess; it is required to resum this large logarithmic term in order to make perturbative 
expansion valid. As for the hi{^'^ / M"^^^^) term in Vl^'^\s)., we emphasize that the coefficient 
of this logarithm turns out to be much smaller than the other terms for a broad range of 
the value r. Especially in the large messenger mass limit this term vanishes. Therefore the 
resummation of the large logarithm is expected to be negligible in 

In order to investigate characteristics of the self-energy functions we consider their 
behavior in a large messenger mass limit (x — )• and Mmcss ^ ^/Fx)- In this limit, they 
become 



f^^'^s) = (j) AAf^ess2C7(/2)C2(G)(^9-41n(-J)) 



(5.10) 
(5.11) 
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(5.12) 



Eq. (5.10) is the well known result used in various literatures for the gluino mass at one-loop 
order. As for il^^) tj^e large logarithm term does not appear in Eq. (5.11) as mentioned 
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above. It should be noted that the factor 9 in Eq. (5.11) is much greater than 2 which are 
shown in Eq. (12) of ref. [15]. It turns out that the difference between the factor 9 and 
2 originates from different treatments of IR behavior. This large factor, as shown in the 
next section, leads to a significant enhancement to the gluino pole mass at two-loop order 
compared with the result in ref. [15]. We retain an IR behavior in the self-energy functions 
at two-loop order while the authors in ref. [15] get rid of an IR divergence. Since the IR 
behavior persists in the self-energy functions O^^^ and H^^^ it is appropriate to keep an IR 
dependence on the gluino pole mass at two-loop order. In this regards, our method for 
the gluino pole mass at two-loop order is more consistent with perturbation theory rather 
than that of ref. [15]. 

6 Numerical Analysis 

In this section we perform a numerical analysis of the gluino pole mass using the three 
self-energy functions in Section 5. In order to illustrate the numerical significance of the 
NLO correction to the gluino pole mass we compare the NLO pole mass with the LO pole 
mass. The strong coupling constant within the standard model is given as as{mz) = 0.118, 
and the Z-boson mass is given as mz = 91.187 GeV, and the top quark mass is given by 
nit = 173.1 GeV. The top quark mass is required for the as running. The final DR values 
of the gauge coupling Og and of the fermion masses can be converted into the MS ones at 
one-loop order ^ using 



Figure 20 shows that the renormalization-scale dependence of the LO and NLO pole 
mass of the gluino, for A = 150 TeV, A/mcss = 200 TeV and A'^mcss = 1- The LO contribu- 
tion to the pole mass has no explicit dependence on the renormalization scale ^t. But the 
introduction of the running gauge coupling brings about the renormalization-scale depen- 
dence on the LO pole mass. Thus the LO pole mass decreases as the renormalization scale 
increases. The scale dependence of gluino pole mass is alleviated at the NLO. 

Figure 21 compares the LO and NLO pole mass of the gluino as a function of the 
messenger mass Mmess, for /i = 1 TeV, A = 150 TeV and N^^ss = 1- For a fixed visible 
supersymmetry breaking scale, both the LO and NLO pole masses are saturated and the 
NLO correction barely changes as the messenger mass scale increases. 

Figure 22 shows both the LO and NLO pole masses increase as the visible supersym- 
metry breaking scale increases for a fixed messenger mass scale. The NLO correction to 
pole mass also monotonically increases as ]\d^Q^^ increases for fixed values of /i and A. All 
in all, the NLO correction to pole mass is roughly 20% of the LO pole masses of the gluino 
among the three plots. 

^The relations between DR and MS schemes at two- and three- loop order are given in ref. [28] and [29]. 
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200 



Figure 20. The renormalization scale-dependence of the LO (dashed hne) and NLO (sohd hne) 
pole masses of the gluino for A = 150 TeV, Mmoss = 200 TeV and A^mcss ~ 1. 



fi = 1 TeV, A = 150 TeV, A^mess = 1 




200 300 
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Figure 21. The messenger scale-dependence of the LO (dashed line) and NLO (solid line) pole 
masses of the gluino for /i = 1 TeV, A = 150 TeV and A^moss = 1. 



-17- 



/i = 1 TeV, M^ess = 200 TeV, iV^ness = 1 




LO 
NLO 



40 60 80 100 120 140 160 180 

A (TeV) 

Figure 22. The SUSY breaking scale-dependence of the LO (dashed hne) and NLO (sohd hne) 
pole mass for = 1 TeV, Af^css 200 TeV and V„icss = 1- 



We consider benchmark planes, lines and points of the MGM in order to interpret the 
experimental results in terms of possible manifestations of SUSY. The authors in ref. [30] 
proposed the benchmark planes, lines and points of the mGMSB model, producing the 
spectra at specific benchmark points that illustrate different possible experimental signa- 
tures. We adopt their benchmark points and calculate the ratio of the NLO pole mass 
correction to the LO pole mass. Table 1 shows the benchmark points along the bench- 
mark line, mGMSBl which is defined by A'^mcss = 3,tan/3 = 15,/i > 0, A = Mmcss/2 with 
AMmess = 10 TeV. Table 2 lists the benchmark point along the benchmark line, mGMSB2.1 
which is defined by Vmess = l,tan/3 = 15,;U > 0,A = 0.9Mmcss with AMmess = 10 TeV. 
We see that the ratio of the NLO pole mass correction to the LO pole mass for the gluino 
can reach 32% in the mGMSBl. The large ratio compared with that of the mGMSB2.1 is 
mainly due to -Vmess ™ O^^^H^^). 

All the numerical analyses indicate that the NLO correction is large enough to reach 
20% of the LO pole mass or even more. There are three ingredients for the large correction: 
(i) the relative strength of the SU (3)c gauge coupling is large, (ii) the color representation 
of the gluino is octet, (iii) and the number of the Feynman diagrams associated with the 
messenger fields are large. The first two reasons are same with that of gravity mediation 
mechanism (i.e. CMSSM) while the last is characteristic of gauge mediation mechanism. 
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Points 


AfrnoQc /TeV 


LO 

CI 


^NLO 








c/UO 


1 956 


39 


mGMSBl.2 


80 


1089 


1436 


0.32 


mGMSBl.3 


90 


1225 


1616 


0.32 


mGMSBl.4 


100 


1361 


1796 


0.32 


mGMSBl.5 


110 


1497 


1977 


0.32 


mGMSBl.6 


120 


1633 


2157 


0.32 



Table 1. Line mGMSBl: A^^ess = 3,tan/3 = 15,// > 0, AM„icss = 10 TcV (masses in GeV). 



Points 


-^^mcss /TeV 




^NLO 




mGMSB2.1.1 


80 


719 


917 


0.275 


mGMSB2.1.2 


90 


809 


1031 


0.275 


mGMSB2.1.3 


100 


898 


1145 


0.275 


mGMSB2.1.4 


110 


988 


1259 


0.274 


mGMSB2.1.5 


120 


1078 


1373 


0.274 


mGMSB2.1.6 


130 


1168 


1488 


0.274 



Table 2. Line mGMSB2.1: TV^css = l,tan/3 = 15,// > 0,A = 0.9Mmoss with AM,noss = 10 TcV 
(masses in GeV). 

7 Conclusion and Outlook 

We have presented the self-energy functions for the gluino of the minimal gauge mediation 
at two-loop order and studied the radiative corrections on the gluino pole mass. The one- 
loop pole mass is the leading order while the two-loop correction is the next-to-leading 
order. The next-to-leading order correction shifts the leading order pole mass by roughly 
20% or even more. This shift is much larger than the expected accuracy of the mass 
determination at the LHC, and should be reckoned with for precision studies on the SUSY 
breaking parameters. 

Not only the gluino mass but also the squark masses are crucial to study the phe- 
nomenology at the LHC. The squark masses also involve the SUSY QCD so that its con- 
tribution from higher order radiative corrections are expected to be large. The numerical 
significance of the next-to-leading corrections to the squark masses deserves a detailed 
investigation which we leave for future study. 

The next-to-lightest-superstmmetric particle (NLSP) in gauge mediation is mostly 
either neutralino or stau depending on the specific regions in parameter space. Using 
the full expressions of the self-energy functions one can evaluate the significance of the 
radiative corrections for the NLSP mass and refine the spectra at the benchmark points of 
the minimal gauge mediation model. 

In this paper we have focused only on the minimal gauge mediation. Other gauge 
mediation models tend to retain the feature that the higher order radiative corrections to 
the gluino pole mass are substantial. Thus when one quantitatively studies complicated 
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SUSY-breaking models associated with gauge mediation one must pay attention to the 
next-to-leading order or higher order radiative corrections to the gluino mass in addition 
to the leading order result. 
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A Master Integrals 

Reduction of Feynman integrals to master integrals is performed in the Euclidean space 
with d = 4 — 2e dimensions. Our convention for the master integral with n-propagators 
{'Pi,V2, .■.,Vn), and /i- loops is 



,7£;/2N/l(4-d) 



1 



d/2j p^'p^...T?^ 



(A.l) 



where Ij are loop momenta. In order to simplify the form of master integrals, several 
parameters are defined as follows, 



X 



y 



iO, r 



m 



2 ' 



m 



2 ' 



^X + i- y/x 

VxT4: + ^/x ' 



m,2 

Pl,2 



1 _^ y^xjx + 4) 

2 2^ ' 

r + x-1 ^ l_^j,2 ^2r{x-l) + {x + 1)2 



2x 
r — r + X 
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lb — y/^2 _ 2r(f — x) + [f -\- x)2 



(A.2) 



where the subscript 1(2) is associated with a plus(minus) sign. 

The master integrals shown in figure 19 are given in terms of the e-expansion as follows, 
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The results of the multiple inverse binomial sums in ref. [31] are used for 132, /as and 134. 

The higher transcendental functions Jij and J^^^^^ are, respectively, the zeroth- and first- 
order terms in the e-expansion of the master integrals. The functions with poly logarithms 
up to second order are given by, 
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We stress that Eqs. (A. 15), (A. 16), (A. 17) and (A. 18) are valid for all values of x and r. 
The constraints on the allowed values of x and r in Eqs. (A.19) and (A.20) are valid for our 
calculation because mi, the messenger fermion mass, is much larger than the gaugino mass. 
We check that Eq. (A. 15) gives a consistent result of ref. [32]. The functions J35(x,r) and 
J37(x, r) are compatible with the result of ref. [33] which uses differential equation method. 
The higher transcendental functions with polylogarithm up to third order are given 

by, 
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As for Jsj (i = 1, 2, 3, 4) we need transcendental functions higher than polylogarithm of 
order 3. Their analytic expressions are beyond the scope of this paper. Instead they can be 
given by a definite integral so that numerical evaluation of them is easily performed. To this 
end, we use the differential equation method in ref. [34]. For example, differentiating /51 
by is represented by a linear combination of other master integrals in a diagrammatic 
way. By solving the differential equation, one can express the function J51 as a definite 



d 
dm? 






integral. 

After all, the functions J5i{i = 1,2,3,4) are given by definite integrals of functions as 
follows. 
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where J7ij are a set of known master integral functions, and ji(x)'s are the master integral 
functions in the limit that all the masses are equal. One can find the expressions for Ji(x) 
in ref. [35, 36]: 

Ji(x) = 21n2(y)ln(l-y) + 6C(3)-6Li3(y) + 61n(y)Li2(?/), (A.25) 
J2M^) = 6C(3) + 12Li3(y) + 24Li3(-y) - 8 ln(y)(Li2(2/) + 2Li2(-y)) 

-21n2(y)(ln(l - y) + 21n(l + y)) . (A.26) 

The coefficient functions Cij{x,r), the factor functions fi{x,r) and the corresponding set 
of master integrals Jij{x,r) are given by, 
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+ Li2 (1 - f) (r - X - 3 - f2(r + x + 1) + 2f(x + 2)) 



1 



+ - (26r2 - r(39x + 82) + x(9x + 25) + 4) + 

r 13 
+ - (26r2 + r(15x + 56) - x(llx + 35) - 82) + ^^^^'^ + x + 1) 



+ T^T ( + X + 1) + f ((r + xf 



12 



1) 



x(x + 4)(f — 1) 
2(r - 1) 

-r (r + : 

x(x + 4)(f — 1) 



r — 1 



+ (r + l)(r + X - 2)^ + ln(r) (r^ + 2r(x + 3) + (x - 3)x - 19) 



.-2/ 1 ^ ^(^ + ~ 1) 



+ 6P(r + x + l) + 



r(x + 2) - x(x + 4) + 2 



r — 1 

In (f) (7^ (r^ + 2r(x + 2) + (x - 4)x - 17) + 6f^(r + x + 1) 



+ In^(r) ( 2f^(r + x + 1) + -f(3r - 2x - 13) + 



1 



x(x + 4)(f — 1) 
2(r - 1) 



- ^(r + l)(r + X - 2)^ + In^ (f) Qf (-r^ + r(5 - 2x) - x(x + 5) - lO) 
3 \ 

+ — r^(r + X + 1) j + 41n(r) In (f) (f(— r + x + 3) — f^(r + x + 1)) 
+ 1±1 iny I f (-r(x - 2) - x2 - 2) + "'^'^ t^^^f ~ - 2r(x + 1) - x + 2 



+ ln2y(f('^i!_il + 2r + ^ + l) + 



r — 1 

(x + 2)(x + 4)(f - 1) 



2(r - 1) 



2x 



(x^ + 2x + 4) 



(x + 2)(x + 6) 2 



+ 
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+ - 1) ln(7?2 X + 1) + (r/2 - 1) ln(7?i x + 1) j i^im - 1) ln(7?2 x + 1) 
+ {fi2 — 1) 111(771 X + 1) ) ( x{—r + X + 3) — xf{r + x + 1) 



+ \ - 1) ln(r/2 X + 1) + (r/2 - 1) ln(r/i x + 1) j ^ - - 7/2) 2; fin (f) 

H ^ + f(2r(r + x - 1) + 6x) + 2 i-rix + 1) + x^ + x + 1) 

r — 1 



+ (^(771 - 1) ln(f/2 X + 1) + (??2 - 1) ln(^i X + 1) 

X 1^ ln(r) {r (r^ + (r + 2)x - l) + +^^)^'^^- ^) _ 2r{x + 1) - x + 2^ 
- 2x2 _ 2 + (A.30) 



- 1 

/4(x,r) = (r - l)~^x~2(r7i - 772)"^ , 

JAj{x,r) = |j^2^(x,r), J35(x,r), J37(x,r), J37^^,^^, J4i(x,r) 

C4i(x,r) = -C45(x,r) = r - x - 1 , 
C42(x,r) = -4, 
C43(x,r) = -r - X + 3, 
C44(x, r) = 3r — X — 1 , 

TT^ 37r 15 

C4o(x, r) = 2(1 - r)Li2(l -r) + — (r + x+l) + — + — - ln(r)(llr + x + 1) 

D 2 2 

+ - ln2(r)(7r - x - 1) + ln(y)(r - x - l)(ln(r) + 1) 

2 1 - y 

+ \T?{y){^—^ + !1±^±1 ] + /"(^^ _ l)ln(?72X + 1) + (7?2 - l)ln(7?ix + 1) 



X {2{r - X - 1) - 2x ^1^^^ ln(y)^ , (A.31) 
where the new variables are defined by 

y' = y\x^x/r^ fll = 'ni\r-^f , m=r]2\r^f- (A.32) 

For the calculation of the gluino pole mass, setting x — )• leads to a good approximation 
of the self-energy function Q^'^K In that limit, we can do the above integrations so that the 
analytic expressions of J^i are given as follows: 

. X (r + 1) , , 4 (3-r) , , , In'^(r) 
J5i(0,r) = -f^-^Li2(l-r) - + ) T^ln(r) 



(r-l)2 ' r-1 (r-l)2 2(r - 1) 

1 ln(r) 
r — 1 (r — 1)^ 



+ ln(x)[:— , (A.33) 



J52(0,r) = J53(0,r) = ^Li2(l-r) + --— 2ln^(r), (A.34) 

r — 1 (r — Ij^ 
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f — 1 r — 1 r — 1 

vr^ rln(r)ln(f) ln(l — f) ln(f) 



+ 1 Iw, ' + 



6(r-l) (r-l)(f-l) f-1 



Li2(J^) +ln(^)ln(r-f) + i^y (A.35) 



It should be noted that there remains a singularity of ln(x) in J^i{x,r) as x — t- 0. We 
comment that the singularity can be handled using the asymptotic expansion technique [37, 
38]. 

B Explicit formula for OS^^s) 

Using the master integrals in Appendix A, we can write down the full expression of the 
self-energy function fi^^^ as follows, 

2 



4j(i)(a;,^)( : ^—^ + 



r-x + 1 ((x + l)2+r(x-5)) 



'22 ' 



+ 4J35( -, - ) I r-x + 2 



{rji - ''I2) "^rx 
{x - 2)(x + 1)2 + r (x2 - 9x + 2) 



r' r J \ ' ' ivi ~ '?2) '^x'^ 

X 1' 



+ 8J36(x, r) + 4J37(x, r) + 4J37 



r r 



2J37 -,- 2r- 3x + 4 + 



2 ((x - l)(x + 1)2 + r (x^ - 7x + 1)) 



r ' r y \ (7/1 — r/2) ^2;2 



^ /X r 
+ 2J37 



(x^ + 2x2 + 3x + r ^^2 _ 73, _ 2) + 2) 



f ' rJ {r]i — 772) 2x 

, , /r-x+1 (x + 1)2 + r(x - 5) 
+ 4J4i(x,r) 



(771 — 772) 2rx 

^ /xl\/ 3r + x + l\ / 3r + x + l 

+ 2J41 -,- h' + 3- 7 +2J42 (x,r,f) l-r + 



r J \ ivi ~ ^2) 2x/ ' ' V (^1 ~ ^2) 2x 

- 4(r + X + 1) J52(2;, r) + 2(r - 1)( J54(x, r, f) - J53(x, r)) 

2Li2(l-r)/, x3 + (r - 4)x2 - (17r + 13)x - 8 , .no 
+ 2 r + 4 X + 3 + ^ 1 ^2 ^1 - ^2 a:2 

- 2Li (1 - f=) f r + ^lliljLll^i:!^^ + + + ^) (^^ + 3x + r(x - 4) + 4) 

^ V - ^2) 2x y 6 V (^1 - m) ^2;2 

2 [rji - 772) ^2;2 15x2 + 88x - 144 

r 2r 
(lla;2 + 80x - 104) (a; + 1)2 + r (llx=^ + 7x2 _ ggQ^ ^04) 

2 (^1 — ^2) 2rx2 

/ r (6x2 _ 232; _ 12) + 3 (2x3 _^ 7^2 _^ 9^ _^ 4\ \ 

+ 2 In r 2r - 12x - 19 + ^ '- 

\ (7/1- 7/2) ^x^ J 
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121 (~) ( '^^ + ^^ + 22;^ + 3a; + r(x2-7x-2) +2 
6 {vi ~ V2) ^2;^ 

Ax^ + Ix^ + 12x + r (4x2 - 31x - 9) + 9 

(x^ - 1) + r (3x2 _ 27a; + le) 



- 2 ln(r) In (f) |^4r + 
+ ln2(r)(2(r + l) + 



(??1 - f?2) ^X^ 

3 (x^ + 2x2 + 3x + r (x^ - 73; - 2) + 2 

(^1 - ^2) 2x2 

x-2 (x + l)2 + r(x-5) 



+ ln2 (f ) \ 2r + 
4(l + y) 

+ n — ^^^yi \ ( \2 

1 + y' . /2(3r + x + l) 



1 + y /2(3r 

+ — 1 My -, — 



ri2) 2x 



4r 



,'2(3x2 + 5rx + 9x + 2r + 6) 8(x + l)\ , 2< ,^( 4 

+ ln2(y) ^ 1 ^T—^ - + ln2(y ) U + - + 22 

rx I \ X 



+ r( 10- - 1 + 

X 



(^1 — '72) 2rx2 
8\ (3x-14)(x + l)2 + r(3x2-41x + 22) 



a;2 (??! - 7?2) ' 



4x (r/i - 772) ■ 



+ 4(^(771 - 1) ln(7/2 X + 1) + (?/2 - 1) ln(r/i x + 1) 



(7?i - 1) ln(r/2 X + 1) + (r/2 - 1) ln(r/i x + 1) 



/ 4(llr-4x + l) 
I r 



+ 



16 ((x + 1)2 + r(x - 3)) _^ / _^ 4 (x2 + 2>x-2r + 2)\ 
im - V2) 2rx V {rji - 772) 2x2 J 



ln(f)( ^ii!-4L_2r + 4) 1 +41n'' ^ 



(??1 - %)2x 



ln(r)( / + " + ! -1 



^mess/ V V(^l -^?2)2x 



+ ( (??i - 1) ln(7?2 X + 1) + (r/2 - 1) ln(??i x + 1) ) ( J^^"^^^^ - 2 



(r -H- f) 



+ (^^J A(r - l)-iiVn,essC(/?)C2(G) 

- 3 J^^) (x, r) - - J35(x, r) + J,-r{x, r) (2 + - - ^ 
X \ X r — 1 

+ J f ^ 1 \ / 8x (771 — 772) ^ 8x2 _|_ _(_ _ 8^ _|_ 8 
yr'ryy r — 1 (r — l)x 

+ 2J4i(x,r) + J42(x,r,f) + J5i(x,r)(2-2r + x) +4J51 ( -,-) ( 1 - - 

\r r J \ r 

+ 2J52(x, r)(r + X + 1) + (r - l)(J53(x,r) - J54(x, r,f)) + 8Li2(l - r) 



■ 1 — r 
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+ 



+ 



+ 




(]iP^))-(^"^' 



(B.l) 
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